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ABSTRACT 

We present 75^' x 75'' size maps of M82 at 6.4 /im, 6.6 /im, 7.7 /im, 31.5 /im, and 37.1 /im with a 
resolution of ^4" that we have obtained with the mid-IR camera FORCAST on SOFIA. We find strong 
emission from the inner 60'' (~lkpc) along the major axis, with the main peak 5" west-southwest of 
the nucleus and a secondary peak 4" east-northeast of the nucleus. The detailed morphology of the 
emission differs among the bands, which is likely due to different dust components dominating the 
continuum emission at short mid-IR wavelengths and long mid-IR wavelengths. We include Spitzer- 
IRS and Herschel/FACS 70 /im data to fit spectral energy distribution templates at both emission 
peaks. The best fitting templates have extinctions of Ay = 18 and Ay = 9 toward the main and 
secondary emission peak and we estimated a color temperature of 68 K at both peaks from the 31 
/im and 37 /im measurement. At the emission peaks the estimated dust masses are on the order of 
10"^ Mq. 
Subject headings: galaxies: ISM — galaxies: starburst — galaxies: individual (M82) — Infrared: ISM 






1. INTRODUCTION 

Over the last few 1 0^ yr a series of starbursts has 
been triggered in M82 (jde Griisll200TI : lMavva et al.ll2006f ) 
due to its interactio n with M81 and NGC 3077 (e.g . 
ADDleton et al. 1981; lYun et al.l 119941 : iSun et al.l 120051 : 
de Mello et al.l I2008D. The mos t recent starbu r st (< 
50 Myr, e.g. iRieke et al.l 119801 : iSatvapal et "all Il997h 
has created a stellar cluster at the center of M82. A 
ring of ionized gas that is enveloped by a molecular 
gas ring surrounds this cluster, and the starburst is 
likely fueled by gas that is funneled toward this re- 
gion by a ^1 kpc (^IQ stellar bar (IJelesco et al. 1991; 
Achtermann fc Lacvlll995l : IWills et"al]l2000l : iGreve et al.l 
20021 ). These gas rings, which appear between about 
10" and 30" from the nucleus, could be gas swept up 
by a central expanding superbubble (jMatsushi ta et al.l 
[2000, 2005). High external gas pressure imposed on the 
molecular clouds probably drives the current starburst 
(JKeto et al. 2005). This view is supported by the de- 
tection of warm molecul ar gas that is likely heated by 
dissipation of turbulence (|Panuzzo et al.|[2QlQ[ ). 

The view to the center of M82 is obscu r ed by 
patchy and high e xtinction (Satvapal et al.l 119951 : 
iForster Schreiber et a l. 2001). The mid-infrared (mid- 
IR) wavelength regime, where the extinction is mod- 
est, is therefore used ex tensively to study the cen- 
ter of M82 (e.g. IRieke et al. 198 0l: iTelesco et al.l [19911: 
IForster Schreiber et al.lMoic iBeirao et al.l 120081 ). Indi- 
vidual star forming clouds that are located at the inner 
edge of the molecular gas ring or within the ionized gas 
ring have been observed in the mid-IR with large ground- 
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based observatories that provide (sub-) arcsecond resolu- 
tion (Lipscy & Plavchan 2004; Gandhi et al. 201 11). Un- 
fortunately, ground-based mid-IR observations are lim- 
ited to a few telluric windows at < 30 jam. 

We have used the mid-IR camera FORCAST on 
SOFIA to study the central region of M82. FORCAST 
provides relatively high spatial resolution and also cov- 
ers the entire 5-40 /im wavelength range. Emission from 
hot, small grains dominates the continuum at short mid- 
IR wavelengths while emission from warm, large grains 
dominates at long mid-IR wavelengths. A total wave- 
length coverage is therefore important to determine the 
detailed contribution of each dust component to the con- 
tinuum emission and to estimate the heating of the dust. 

We adopt a distance of 3.5 Mpc for M82 
(jDalcanton et al.l l2009f ). so that 1" corresponds to 
17.5 pc. 

2. OBSERVATION 

We observed M82 in the 6.35 /im (AA = 0.14/im), 6.61 
/im (AA = 0.24/im), 7.71 /im (AA = 0.47/im), 31.5 /im 
(AA = 5.7/im), and 37. 1 /im (AA = 3.3/im) passband 
filters using FORC AST (lAdams et alllMol : [Herter et al.l 
[20T2bh on SOFIA (jYoung et al.ll2012D . 

FORCAST is a dual-channel camera with a short- 
wavelength channel (5 — 25 /im) and a long- wavelength 
channel (25 — 40 /im), allowing for simultaneous obser- 
vations at two specific wavelengths by selecting band- 
pass filters in filter wheels. A dichroic beam-splitter di- 
rects the same field of view to both the short- and long- 
wavelength cameras. To increase sensitivity, the beam- 
splitter can be either removed or replaced by a mirror 
during an observation. The effective field of view of the 
two 256 square arrays is 3.4' x 3.2^ After distortion cor- 
rection, the pixel size is 0.768''. 

We observed M82 on 2010 December 1 in the two 
long wavelength bands without dichroic at an altitude 
of about 43000 feet. We used the on-chip chop- nod ob- 
serving mode with a chopper throw and nod-distance of 
120" and 90", respectively, and a chop frequency of 2 Hz. 



Nikola et al. 



In this mode the source is alway on the detector array. 
For the 37 /im and 31 /im observations we took seven and 
five integrations, respectively, with an on-source time of 
60 s per integration. The observations in the 37 /im band 
were taken at a zenith angle of about 58°. 

On 2010 December 4, we obtained four integrations of 
M82 in the 6.6 /im and 7.7 /im bands and eight integra- 
tions in the 6.4 /im band from an altitude of about 42000 
feet. All bands were observed without beam-splitter, ex- 
cept for the 7.7 /im band observations, which were ob- 
served simultaneously with the 31 /im band. We ob- 
served again in the on-chip chop- nod mode, but with a 
chopper throw and nod-distance of 90'' and a chop fre- 
quency of 5 Hz. All bands were observed at a zenith 
angle of about 53°. 

Since M82 is bright in all bands, we spatially regis- 
tered the individual chop/nod integrations using a two- 
dimensional Gaussian fit over the entire galaxy before 
CO- adding the individual integrations. These pointing- 
corrections are about 3 pixels and improve the noise and 
the final maps. Due to the on-chip chop- nod observing 
mode the extent of the final images is about 75'' x 75". 
Since pointing accuracy and drifts were on the order of 
several arcseconds, we determined the absolute position 
of the M82 observation by comparing the 7.7 /im FOR- 
CAST map with the 8 /im Spitzer/IRAC map from the 
Spitzer archive. 

The flux densities in all bands were color-corrected. 
The standard stars /3Gem, /3UMi, and /iUMa were used 
for absolute flux density calibration and to estimate the 
beam size. Based on the uncertainty of flat fielding and 
the water vapor burden, as well as the variance of the 
flux density measurements, the uncertainty of the in- 
trinsic flux density of the standard stars, and the fact 
that /3Gem and /3UMi are variable stars, we estimate 
a 3a uncertainty of ^20% for the absolute calibration 
(jHerter et al.ll2012al lbl). Due to turbulent airflow across 
the telescope and pointing instabilities the observations 
are not diffraction limited. For the M82 observations, we 
estimate a beam size of about 4" ±0.5" for all the bands. 

The original maps of the 6.4 /im, 6.6 /im, and 7.7 /im 
continuum had a lower signal-to-noise than the two long- 
wavelength maps and we therefore smoothed them with a 
2 pixel FWHM Gaussian kernel. The final flux densities 
are presented in Table [H 

We include Hers chel/F ACS 70 /im photometric data 
in our spectral energy distribution (SED) analysis below. 
Since processed and published PACS 70 /im photomet- 
ric data were unavailable, we used level 2.5 data product 
from the Herschel archive. The observations were car- 
ried out in scan-map mode at medium speed (20" s~^) 
resulting in a smeared beam of about 5.46" x 5.76" (PACS 
Observer Manual). We extracted background-subtracted 
flux densities within a 6.8" aperture (~ 3.5 PACS pix- 
els) at three positions (see Table [1] for positions). The 
flux densities are 115 ± 11 Jy (main peak), 98 ± 10 Jy 
(secondary peak), and 116 ± 11 Jy (ridge). Within a ra- 
dius of 35", similar to the entire region observed with 
FORCAST, the 70 /im flux density is 1644 ± 41 Jy. 



3. RESULTS 
3.1. Morphology 



Figure [T] shows the distribution of the mid-IR emis- 
sion in each of the observed bands. The 6.4 /im, 6.6 
/im, and 7.7 /im emission is slightly more extended than 
the 31 /im and 37 /im emission. Also, the contrast 
in the morphology is higher in the 31 /im and 37 /im 
emission than in the short-wavelength bands, with the 
peaks being more pronounced. The main peak of all 
mid-IR distributions is 4.5" west-southwest from the dy- 
namical center. Strong emission extends further to the 
west-southwest and "peaks" roughly 9" from the center 
("western ridge"). A secondary peak is visible 4" east- 
northeast of the center (except at 6.6 /im). In the 7.7 /im 
band, the emission at the secondary peak appears as two 
separate components. The positions of the main and sec- 
ondary peaks coincide roughly with the position of the 
ionized ring or the inner edge of the molecular ring. 

The 7.7 /im FORCAST map is in excellent agreement 
with the Spitzer/IRAC 8 /im map (Figure [T]). Also, the 
FORCAST 31 /in i map agrees v ery well with the 30 /im 
map obtained by iTelesco et al.l (jl991) using the NASA 
Infrared Telescope Facility (IRTF), and the flux densi- 
ties at the main and secondary peaks agree within the 
calibration uncertainty. 

In Figure [2] we show a three-color image of M82, where 
we have combined the 6.6 /im (blue), 31 /im (green), and 
37 /im (red) bands. All colors are scaled linearly with 
flux density and start at the 3<j level of the statistical 
background noise in each band. The image shows that 
the 6.6 /im is more extended and that the secondary peak 
and western ridge are more pronounced in the two long- 
wavelength bands. 

Figure [3] shows the emission profiles along a 60" cut 
following the major axis of M82. The length and position 
of the cut is shown in Figure [TJ The profiles are in steps of 
1 pixel and the flux densities are summed over 1x5 pixel 
strips perpendicular to the major axis and normalized 
to the emission peak. They emphasize that in the 31 
/im and 37 /im bands the ratio between the main and 
secondary peak is much stronger than in the three short 
mid-IR bands. In fact, the 6.6 /im continuum only shows 
a gradient in the flux density from the main peak toward 
the position of the secondary peak. The emission from 
the "western ridge" is very pronounced in the 31 /im and 
37 /im bands, less in the 7.7 /im and 6.4 /im maps, and 
barely noticeable in the 6.6 /im map. The 31 /im and 37 
/im emission show a peak at this position that is stronger 
than the secondary peak. 

Figure [3] also shows the flux density ratios along the 
major axis of M82. The profile of the 37/im/6.6/im and 
37/im/7.7/im fiux density ratios basically follows the pro- 
file of the 37 /im flux density. In contrast, the flux density 
ratios of the 7.7/im/6.6/im and 7.7/im/6.4/im bands are 
fairly uniform along the major axis, with larger variations 
only at the edge of the profile, where the signal-to-noise 
becomes small. The profile of the 37/im/31/im flux den- 
sity ratio differs from all others. It is basically flat within 
±10" of the center of M82 and then increases on either 
side, while the flux density ratios of all other bands either 
decrease or stay constant with distance from the center. 
This indicates a higher color temperature of the large 
grains in the inner region than in the outer region. Also, 
the difference of the profiles between the three short mid- 
IR bands and two long mid-IR bands indicate that they 
trace different dust components with spatially distinct 
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excitation conditions. 

3.2. Mid-IR Spectral Energy Distribution 

iSiebenmorgen fc Kriigell (|2QQ 7l)FI have created template 
SEDs for starbursts as a function of the total luminos- 
ity (I/(tot)), optical extinction, contribution to the total 
luminosity by OB stars, and radius of the nuclear star- 
burst, covering the wavelength range 0.03-2000 /im. We 
use the FORCAST and PACS 70 /im data points as well 
as low-resolution Spitzer-1Y{S spec tra between 5.3 /im to 
12.8 /im from region 2 and 3 of iBeirao et al.l (|2008l ) to 
find the SED template that fits the data in the two mid- 
IR emission peaks within a 6.8'' x 6.8'' (~ 115 x 115 pc) 
area, corresponding to ^9 x 9 FORCAST pixels. The 
lower size limit of the area is given by the Spitzer obser- 
vation. Since this area is much smaller than the smallest 
nuclear radius in the SED template grid (350 pc), we 
constrained our search algorithm to that template series 
and allow it to fit an additional scaling factor. Figure H] 
shows the data points and the best-fitting SED template 
at the main and secondary peaks. It includes three data 
points from iTelesco et al.l ()l991f ). multiplied by a fac- 
tor of two, which corresponds to the flux density ratio 
within a 6.8" and 4.5" aperture in a FORCAST band at 
the peaks. The flux densities at 70 /im and 37 /im at the 
main and secondary peaks are very similar and the 70 
/im flux density is thus crucial to constrain the SED fit. 

For the main mid-IR emission peak, the best-fitting 
SED template has a total luminosity of LMp(tot) = 
10^^"^ 1/0, an OB luminosity fraction of 40%, an extinc- 
tion of Ay = 18, and a hydrogen density of n = 5 x 10^ 
cm~^ in the hot spots. The template requires an addi- 
tional scaling factor of 0.041, which is of the order of the 
ratio of observed to the template area, and lowers the 
total luminosity to LMp(tot) = 6.7 x 10^ Lq. The best 
fitting model for the secondary emission peak has a total 
luminosity of Lsp(tot) = 10^^-^ L©, an OB luminosity 
fraction of 40%, an extinction of Ay = 9, a hydrogen 
density of n = 5 x 10^ cm~^, and requires an additional 
scaling factor of 0.11, which reduces the total luminosity 
to Lsp(tot) = 5.7 X 10^ Lq. 

Due to limited resolution of the SED template param- 
eters we estimate the derived extinction to be within 
a factor of two. B oth values are lower th a n the ex- 
tinction reported by iForster Schreiber et al.l (|2001f ) for 
the entire central region {Ay = 52 ± 17). This dis- 
crepancy could be due to the patchy extinction in the 
central region ([SatyaDal et al.l [19951 : iLipscy fc PlavchanI 
[200llGandhi et al.ll201lD " 

3.2.1. Extinction 

The extinction laws between 3 and 10 /im toward the 
cent er of M82 and the Galactic center (GC) are simi- 
lar (jForster Schreiber et al.ll2001|). and for the GC it is 
best modeled with Ry = 5.5 ()Draiii3l2003f ). To estimate 
the extinction and emissivity at the FORCAST bands 
we therefore use the GC extinction law as modeled by 
iLi fc Drainj (2001) and Weingartner fc Draine (2001)0 
with Ry = 5.5, a gas-to-dust mass ratio of 105, a dust 
grain mass per hydrogen atom of rrid = 2.2 x 10~^^ g/H, 

6 www.eso.org/~rsiebenm/sb_models/ 

^ www.astro.princeton.edu/~draine/dust/dustmix.html 



N{R)/Ay = 1.37 X 10^^ cm-2 mag-^ (with N{R) = 
A^(HI) + Ar(H2)), and an emissivity index of /3 = 1.79 for 
the 31 /im and 37 /im bands. The mass absorption coef- 
ficient and albedo at 31.62 and 37.15 /im are hi = 385.2 
cm^g-i and 0.0007, and tz = 288.8 cm^g-i and 0.0006, 
respectively. Using the relation K,ahs{^) = 0-4 In 10 (1 — 
albedo) x [A{X) / Nu] / m^ust cm^g-i (|Li fc Drainell2Q01f ) 
leads to A{37)/Ay = 9.46 x 10"^ and A{31)/Ay = 1.26 x 
10-^ resulting in r(37/im) = 0.16 and r(31/im) = 0.21 
at the main emission peak, respectively. At the sec- 
ondary emission peak the opacities are r(37/im) = 0.08 
and r(31/im) = 0.11. 

3.2.2. Dust Mass 

We use three different methods to estimate the dust 
masses at the main and secondary emission peaks within 
6.8'' X 6.8''. 

Applying the extinction laws from the previous section 
we estimate the gas mass using Mg = /i x mn x A^(H) x 
Area, where /i = 1.4 is the mean atomic mass per hy- 
drogen, mu is the mass of atomic hydrogen, and A/'(H) 
is the total hydrogen column density. We obtain a total 
gas mass of Mg(MP) = 3.8 x 10^ M© at the main peak 
and Mg(SP) = 1.9 x 10^ M© at the secondary peak. Us- 
ing a gas-to-dust ratio of 105 this gives dust masses of 
Md(MP) = 3.6 X 10"^ Mq and Md(SP) = 1.8 x 10"^ M© at 
the main and secondary peak, respectively. Due to the 
uncertainty of the extinction, the dust masses can vary 
by a factor of two. 

The dust mass can also be estimated by M^ = l/^^abs x 
Fj^{X) X D'^ / B{\^T)^ where Ai:abs is the mass absorption 
coefficient, Fj^ is the ffux density, D is the distance to 
M82, and B{\^T) is the Planck function at wavelength 
A and temperature T. Using the previously estimated 
opacities for 31 /im and 37 /im and a modified blackbody 
function we estimate a color temperature of 68 ± 10 K 
at both emission peaks. The resulting dust masses are 
1.15 X 10"^ M© and 7.9 x 10^ M© at the main and 
secondary peaks, respectively. The uncertainty of the 
masse s is within a fact or of four. 

San ders et al.l (|1991l ) estimate the dust mass using the 
relation M^ = (L(FIR)/10^ L©) x (40K/Td)^ 10"^ M©, 
where T^ is the dust temperature and L(FIR) is the 
far-IR luminosity in the range between 40 and 500 
/im. From the best fitting model SED we determine 
a Lmp(FIR) = 3.08 x lO^L© at the main peak and 
Lsp(FIR) = 2.52 X lO^L© at the secondary peak. Us- 
ing the color temperatures derived above we obtain dust 
masses of Md(MP) = 2.2 x lO'^M© and Md(SP) = 
1.8 X lO^M© at the main and secondary peak. These 
dust masses are good to a factor of three. 

High resolution CO measurements with a 4.2" aperture 
suggest molecular hydrogen column densities of A/'(H2) ~ 

(4 10)xlO^^ cm~^ at our main peak and A^(H2) ^ 4x 

10^^ cm~^ at our secondary peak (jWeifi et al.ll2001l ). So, 
the enclosed gas masses, including He and other heavy 
elements, are (1.7 4.2) x 10^ M©. Assuming a gas-to- 
dust ratio of 105 this yields dust masses between (1.6 — 
—4) X 10^ M©, in agreement with our calculations above. 

4. SUMMARY 

We have presented the first results of mid-IR observa- 
tions of M82 obtained with FORCAST on SOFIA. M82 
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was observed in the 6.4 /im, 6.6 /im, 7.7 /im, 31.5 /im, 
and 37.1 /im bands. The observations cover a 75" x 75" 
region. Ah bands show a strong peak located 4.5'' west- 
southwest of the kinematic center of M82. A secondary 
peak 4" east-northeast of the nucleus is seen in the 6.4 
/im, 7.7 /im, 31.5 /im, and 37.1 /im bands, but not in 
the 6.6 /im continuum. The profiles of the flux density 
ratios over 60'' along the major axis indicates that the 
emission at the three short mid-IR bands is dominated 
by a different dust component than the emission at the 
two long mid-IR bands. We fitted SED templates to 
the FORCAST data combined with PACS 70 /im and 
low-resolution Spitzer-IRS spectra and estimated extinc- 
tions of Ay = 18 and Ay = 9 toward the main and 
secondary peaks and a dust color temperature of 68 K 
in both peaks. The dust masses at the locations of the 
emission peaks are on the order 10^ Mq. 
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sion Laboratory, California Institute of Technology un- 
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Table 1 

M82 FORCAST Flux Densities 



Band 


^peak 


C b.c 
'-'mp 


o b,d 
'-'sp 


C b.e 
Owr 


q f,C 


o f,d 
'-'sp 


S'(Total Map)g 


(/im) 


(Jy pixel" -'^) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


6.4 


0.112 ±0.007 


3.42 ±0.23 


2.96 ±0.20 


2.75 ±0.18 


6.93 ±0.46 


6.02 ±0.40 


68 ±5 


6.6 


0.047 ±0.003 


1.40 ±0.09 


1.21 ±0.08 


1.17 ±0.08 


2.90 ±0.19 


2.52 ±0.17 


32 ±2 


7.7 


0.141 ±0.009 


4.12 ±0.28 


3.41 ±0.23 


3.20 ±0.21 


8.16 ±0.54 


6.94 ±0.46 


75 ±5 


31.5 


1.86 ±0.12 


57.1 ±3.8 


40.6 ±2.7 


46.5 ±3.1 


110.0 ±7.3 


79.7 ±5.3 


676 ± 45 


37.1 


2.42 ±0.16 


74.5 ±5.0 


51.5 ±3.4 


61.4 ±4.1 


143.9 ±9.6 


102.2 ±6.8 


891 ± 59 



^ Pixel size: 0.768^^ 

^ Within 6x6 pixels, corresponding to a 4.6^^ x 4.6^^ region. 

^ Main peak: 09^55^51.28^ ±69°40'45.5'' 

^ Secondary peak: 09^55^52. 68^ ±69°40'48.5'' 

^ Western Ridge: 09^55^ 50.47^ ±69°40'43.9'' 

^ Within 9x9 pixels, corresponding to a 6.8^^ x 6.8^' region. 

s Within 50^^ x 75^^ (65 x 98 pixel) region around center of M82 
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Figure 1. Maps of M82 in the FORCAST bands in units of Jy pixel-^ at 6.4 /im (top left), 6.6 /im (mid left), 7.7 /im (bottom left), 31.5 
/im (top right), and 37.1 /im (mid right). Bottom right: IRAC band 4 (8 /im) map (in MJy sr"-*^) overplotted with FORCAST 7.7 fim 
contours . The color scale is linear and starts at the 3cr level of the statistical background noise (0.009 Jy at 6.4 and 6.6 /im, 0.018 Jy at 
7.7 /im, 0.042 Jy at 31.5 /im, and 0.051 Jy at 37.1 /im). The dashed line in the 6.6 /im map (mid left) indicates the position of the profiles 
shown in Figure [3] 
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Figure 2. Three color image of M82, with the FORCAST 6.6 /im as red, 31.5 //m as green, and 37.1 /im as blue. All bands are linearly- 
scaled, starting from 3a of the statistical background noise, before combined. 
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Figure 3. Flux densities (left) and flux density ratios (right), normalized to the value at the main peak position, along the major axis of 
M82. The reference position is the main peak, distances are in arcsec, and positive distance is toward the northeast. Flux densities are 
summed over 1x5 pixels perpendicular to the major axis. 
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Figure 4. Mid-IR SEP of the main peak (left) and secondary peak (right). The sohd (black) line is the low-resolution Spitzer/IKS 
spectrum (jBeirao et al.ll2008l V filled triangles (red) are the FORCAST observations, filled diamond (blue) is the Herschel/PACS 70 /im 
observation, filled squares (green) are IRTF observations (Telesco et al. 1991), multiplied by a factor of two (see the text), dashed line 
(blue) is the Siebenmorgen & Kruegel SED model ^i ebenmorgen Sz KrUge l ^2007). The error bars are smaller than the symbols. 



